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We investigate the magnetic field distribution inside a Type II superconductor which has point 
group symmetry O such as Li2Pt3B. The absence of inversion symmetry as a departure from perfect 
cubic group Oh causes a magnetization collinear with the phase gradient associated with the order 
parameter, and a component of current collinear with the local magnetic field. In the vortex state, 
we predict, by solving the Maxwell equation, a local magnetic field transverse to the vortex lines. 
The probability distribution of this transverse field is also obtained for the prospective muon-spin- 
rotation measurements. 
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I. INTRODUCTION 

Noncentrosymmetric heavy-fermion superconductors 
have recently triggered much theoretical and experimen- 
tal attention. The family of Li2B(Pdi_2;Pta;)3 as well as 

\ the compounds CePtsSi and CeRhSis are among the ma- 
jor examples. i^'^i^'^'^'i'^'^ The lack of inversion symme- 

I try of background crystal causes an unconventional pair- 
ing in the superconducting state and a number of conse- 
quences are predicted] It is commonly believed 
that in the normal state the Fermi surface is split by the 
spin-orbital coupling as a result of the broken inversion 

\ symmetry, and the pairing wavefunction is in a mixture of 
singlet and triplet pseudospin state.™ Suggested by the 
low temperature measurements of penetration depths, 
the family of Li2B(Pdi_3;Pt2;)3 with point group sym- 
metry O demonstrates a trend from a conventional gap 
to a nodal gap caused by mixing with a stronger triplet 
pairing as the composition x varies from to li^ The con- 
clusion that the Li2Pd3B compound is a fully gapped su- 
perconductor is consistent with other measurements such 
as NMRjii specific heatfi^ and muon spin rotation (/iSR) 
measurements ji^iii The nodal gap structure of Li2Pt3B 
hinted by the penetration depth measurement, however, 
is controversial since the specific heat and /iSR measure- 
ments suggest conventional pairing. ^^'^^ 

Microscopically, the lack of inversion symmetry allows 
a parity-breaking spin-orbital interaction which leads to 

[ a magnetization induced by the phase gradient asso- 
ciated with the superconducting order parameter, and 

' vice versaJ^ For an ideally two-dimensional system, the 
phase gradient induced by the constant magnetic field, 
or the so-called helical phase, might be measured using 
a Josephson junction between a conventional and a non- 
centrosymmetric superconductor However, this scheme 
is not feasible when considering the decay of magnetic 
field in a real three-dimensional bulk superconductor. On 
the other hand, a unique and macroscopic distribution of 
magnetic field is generated inside the noncentrosymmet- 
ric superconductors due to the induced phase graident 
and magnetizationii2ii2£i2ii^ Therefore, the measurement 
of the local field distribution, such as the /zSR^^ can be 



used to detect the unconventional pairing. Given the fact 
that these noncentrosymmetric compounds are strong 
Type II superconductors f^i^ the effects of lacking inver- 
sion symmetry on the vortex state are of importance. In 
this paper we focus on the noncentrosymmetric supercon- 
ductors with point group symmetry O, such as Li2Pt3B, 
where a distribution of transverse magnetic field is gen- 
erated in addition to the ordinary field parallel with the 
vortex line.^ 

In Secini we first derive the magnetic distribution as- 
sociated with an isolated vortex line. In particular, we 
note a singular bound current defined by the curl of mag- 
netization appears within the vortex core. It appears be- 
cause the magnetization here is proportional to the gradi- 
ent of order parameter phase associated with the vortex 
state. Hence a new boundary condition is necessary for 
the azimuthal field to balance the singular bound current, 
which was missing in our previous study^^S and it leads 
to a significant correction to this field. In Sec lIIIl we ex- 
tend to the lattice and numerically obtain the tranverse 
field Bi^{x, y) as well as the associated probability distri- 
bution p{B±). In Sec II VI we connect the present results 
with the proposed helical phase, and draw the conclusion 
in the end. 



II. ONE ISOLATED VORTEX LINE 

Inside a type II conventional superconductor in the 
presence of external magnetic field, the Abrikosov vortex 
state with a periodic distribution of internal magnetic 
field is developed when the applied field is larger than 
Hci, below which a complete Meissner effect occurs, and 
less than Hc2, beyond which the superconductivity only 
appears on the surface and completely disappear for fur- 
ther increased field. The internal magnetic field inside 
a conventional superconductor is parallel with the exter- 
nal one, whereas a transverse component is developed 
inside this noncentrosymmetric superconductor, as will 
be shown below. Since the internal field for an isolated 
vortex is essential to the vortex lattice, below we first 
solve the single vortex case. 
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Consider a noncentrosymmetric superconductor with 
an isolated vortex line along the z direction. Below the 
cylindrical coordinate {p, (p, z) is employed. We denote 
the order parameter associated with the single vortex as 
\ip\e'^'^ where the phase x=— </>, that is the vortex is said 
to have the winding number of minus one. Besides, the 
magnitude can be viewed as a constant for positions 
outside the vortex core that roughly has a size of coher- 
ence length ^. The internal local magnetic field B is then 
determined by the supercurrent J through the Ampere's 
law, 



V X B = 47rV X M 



47r 



(-e)J, 



(1) 



where the first term on the right represents the bound 
current due to the magnetization M, which is zero for 
conventional superconductors. The electric charge has 
been written as (~e). Here the nonzero M is proportional 
to the phase gradient ji^iii^ namely. 



4^M = 4(A + :f Vx) 



(2) 



where A is the London penetration depth, b has the di- 
mension of length, and the small number (5/ A = k/2 de- 
pends on the strength of the spin-orbital interaction as- 
sociated with the lack of inversion symmetryi^ On the 
other hand, the internal magnetic field also induces a con- 
tribution to the supercurrent, which together with the 
ordinary contribution can be written as. 



-eJ^-(A + -Vx)-^B 



(3) 



where A is the vector potential defining the internal field. 
Note that ?i=l for convenience. M in eq ([2]) and the 
last term in eq. ([3]) exist only if the superconductor is 
non-centrosymmetric. We note here that^^ both 1/A^ 
and (5/A^ are proportional to 1 — 1^(T), where ^(T) is 

the Yosida function, so k is oc (1 — j near and 

increases with decreasing temperature, saturating to a 
constant value as T ^ 0. 

The solutions to Eq. ([1]) go back to the standard singly 
quantized vortex when k is set to be zero. The corre- 
sponding internal field — B^z^ = i^^Kf){p/ X) where 
$0 is a quantum of flux and is the modified 
Bessel function of zeroth order. For nonzero k, an az- 
imuthal internal field B^ emerges. In addition, near the 
core the magnetization contains an azimuthal component 
proportional to 1 / p from the phase gradient term in eq. 
([2|) . This diverging behavior of results in a diverging 
B^ (X which was unfortunately left out in our ear- 
lier worki"^^ To see this in more detail, we first note that, 
due to the Meissner effect, we expect that B, the gauge 
invariant velocity (Vx + if A) hence also M, J all decay 
to zero as p ^ oo. We note that this guarantees the total 



magnetic flux associated with a singly quantized vortex 
remains to be $o, as §{A. + ^Vx) • = over a large 
loop encircling the vortex. Integrating eq ([1]) over a large 
area enclosing the vortex and using Stokes theorem, we 
obtain 



r -.An 
B-dl-A-K (bM-dl^ — (-e) 



J dS 



(4) 



hence automatically the condition of zero total 
current along the vortex line, since B^, M^f, 
zero exponentially at large p. 

However, considering eq ^ for a small loop 
the vortex, we see that for p ^ 0, B^p — AttM^. 
([2]) and the fact that the magnetic flux passin, 
a small loop must necessarily approaches zero 
size of the loop does, we get A-kM^ x 27rp = —4 



transport 
decays to 

encircling 
Using eq 
g through 
when the 
$0 hence 



B^{p ^ = - 



k_%_ 
2X2np 



(5) 



This behavior can also be obtained directly by sub- 
stituting eq and ([3]) into ([T]) and using V x Vx = 

-27r5(2)(-)£^ YVe get 



V X B 



c 



K 

2A 



B-4',5(2)(p)i 



(6) 



The delta function can be taken care of by an appropriate 
boundary condition associated with the azimuthal field 
B^ at the origin, as in eq ([5|). 

Eq. ([5]) can be solved in the manner of perturbation in 
K. For p 0, the delta function can be ignored and the 
first-order correction i?^^'' satisfies the following equation, 
obtained by taking the curl of eq. ^ and using eq. dSj, 



d ^ d 1 
dp p dp A^ 



K d 
X dp 



(7) 



in which the source term on the RHS arises from the extra 
supercurrent induced by the unperturbed internal field 
Bf^ as well as the second term in eq. ([B]) . This equation 
is identical to eq. (37) in our previous worki^ There 
the inhomogeneous solution satisfying B^ — > near the 
origin was obtained and we shall denote it as -^^fix) 
with the argument x defined hy x = p/X. (Thus f{x) 
represents the right hand side of eq (41) in Ref. 22). We 
note that this inhomogeneous solution does not satisfy 
the boundary condition in eq. ([5]), which can be fixed 
by adding an appropriate homogeneous solution to eq. 
([7]) satisfying eq ([5]). Since Ki{x) fulfills eq. ([7]) in the 
absence of the source term and Ki{x) (x 1/x for small x, 
we can easily write down the desired azimuthal field as, 



27rA2 



fix) - -K,{x) 



(8) 
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The azimuthal field still decays exponentially at infinity, 
justifying the previous assumptions. However, in con- 
trast to our previous results,— due to the Ki term the 
azimuthal field for an isolated vortex line now has in- 
creasing magnitude towards the core. 



III. LATTICE OF VORTEX LINES 

Since those heavy-fermion superconductors usually 
have low transition temperatures and large Ginzburg- 
Landau parameter, Hd is practically small and usually 
much less than the external magnetic field, which causes 
the superconductor to be in the vortex lattice state in 
most cases. The components for transverse magnetic 
fields Bx{x, y) and By{x, y) inside the superconductor are 
given by the vector sums of eq. ^ over the vortex lines 
located at each of the lattice vertex. For the family of 
Li2B(Pdi_j;Pt2,)3, the ratio of A/^ is about* 20 indepen- 
dent of the Pt composition. The transverse field depends 
on the lattice geometry, and d/X is the only parameter 
with d being the lattice constant. Figs. [T] and [5] show the 

contours of transverse field = \J B'^ + B'^ in unit of 

Bq — K$o/27rA^ with (i/A=0.5 for a square and trian- 
gular lattices, respectively. Due to the delta function in 
eq. ©, the regions close to the core have larger trans- 
verse fields. The destructive interferences from neighbor- 
ing vortex lines occur at the middle points on each side 
and at the center of each square and triangle, respec- 
tively. 

Statistics on the local internal field B±{x,y) give the 
probability distribution p{B±/Bo). In Fig. [3] and Fig. 
m we obtain the probability distributions p{B±/Bo) for 
various values of d/X in the square and triangular lat- 
tices, respectively. p{B±/Bq) vanishes for B± 0, and 
rises up linearly with B± for small B±. It reaches a 
peak corresponding to the longest contour in Fig. [T] or 
O then decreases with increasing B±. As the value of 
d/X decreases, the vortex lines are getting denser and 
the magnetic fields associated with the peak of p{B±) 
become greater. One should note, however, that there 
is a cutoff magnetic field associated with the vortex core 
for each value of d/X because the minimum length scale 
for the above theory to be applicable here is the coher- 
ence length ^. The colored arrows are used to label these 
cutoffs for various d/X. That is, the values of B± beyond 
these arrows are un-physical and the actual values are 
expected to be much smaller. We note that the cutoff for 
a single vortex is approximately given by ^Ki{^/X) w 10 
times Bq. For low density of vortices, i.e. d/X ^1, the 
cutoff remains the same as that of a single vortex. For 
higher density cases, these cutoffs are modified due to the 
neighboring vortices. 

The distributions p{B±) can in principle be obtained 
from /xSR measurements, operated in the following man- 
ner. Positive muons with spin polarizing along the exter- 
nal field and hence the vortex lines are injected into the 



superconductor in the vortex state. In such a set-up for 
conventional superconductors, the muon spin polariza- 
tion is not affected by the local magnetic fields, as they 
are parallel (assuming that the vortex lines are straight, 
c.f. below). The decay product, the positrons, has an 
anisotropic distribution, being preferentially along the 
initial muon spin. However, for our non-centrosymmetric 
superconductor, the muon spins will precess around the 
local field B±{x,y). This gives rise to an extra decay or 
oscillation of the anisotropic positron distribution. This 
distribution is further a function of external field and 
temperature, with magnetic field magnitude increasing 
with decreasing temperature if below Tc- These features 
may be used to distinguish from other contributions. 
Thermal fluctuation of vortex lines may also give rise 
to horizontal fields, but they would increase instead of 
decrease with temperature. Magnetic fields from nuclear 
moments would have a temperature dependence uncor- 
related with T^,. 

The values of d/X correspond to average external mag- 
netic fields Bav according to the relation -Bau = ~-^ 2jA^ ' 
where 77 is a factor depending on the lattice geome- 
try. For square lattice rj—l, and for triangular lattice 
77=-^. The magnetic field -^j^ is about 25 Gauss given 
A(T 0)=360nm in Li2Pt3B4 If the parameter k ex- 
ceeds 10~^, then the typical fields according to Fig[3]and 
Ulwill be of order 1 gauss, which is supposed to be mea- 
surable in the /iSR measurement. ^"^'^^ In the compound 
Li2Pt3B, however, k may unfortunately be of order 10"'^ 
only, with transverse fields that are hard to detect. Nev- 
ertheless, our calculations apply to any noncentrosym- 
metric superconductor with O symmetry. The transverse 
fields may be easily measurable if k is a few times larger 
and/or A a few times shorter. 



IV. DISCUSSION 

The above paragraphs demonstrate the transverse 
magnetic field in the vortex state as a signature of broken 
inversion symmetry. Here we elaborate on the general 
phase gradient Q = Vx + since this vector is as- 
sociated with the helical phase which is described by an 
order parameter j^'je^'^ ''. In a strictly two-dimensional 
system with open boundary, the phase gradient can be 
obtained by directly setting the current density in eq. ([3]) 
to be zero. Since there is no Meissner effect in such an 
ideal system, Q is determined solely by the external field 
and was argued to be measurable by a Josephson junc- 
tion between a conventional and a noncentrosymmetric 
superconductor,'' 

For our three-dimensional superconductor, Q can be 
obtained from eq. ^ and ([3]) given the internal magnetic 
field B. For positions not exactly at the vortex cores, 

= XkB- X^{V xB) . (9) 
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In a real bulk superconductor, it can be seen that the 
phase gradient is small and not appreciable due to the 
Meissner effect. For example, consider a noncentrosym- 
metric superconductor occupied by a single vortex. The 
unique phase gradient Qz is along the line. Given the 

internal fields si"'' and S^^' obtained previously, Qz is 
vanishingly small at positions p 3> A. Only near the 
vortex Qz is large and given by ^^^j^ for p <C A. We 
note however that, at these distances, Q4, — ^, hence 
\Qct>\ ^ \Qz\- The phase gradient associated with the 
helical phase in this noncentrosymmetric superconduc- 
tors has also been obtained beyond the present modified 
London theoryi^ (see alsoi). 



such as Li2Pt3B by calculating the transverse magnetic 
field in the vortex state. The probability distribution of 
the field is also calculated for the prospective /LtSR exper- 
iments as a signature of the broken inversion symmetry. 
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FIG. 1: Contour plot of transverse magnetic field magnitude 
B± for a square lattice. The unit of B± is Bo = K'l>/27rA^. 
The distance d between neighboring lines to the penetration 
depth A is 0.5. 
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FIG. 2: The same contour plot as Fig. [T]but for a triangular 
lattice with the same value of d/\. 
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FIG. 3: The probability distribution p{B±/Bo) of internal 
transverse fields for square lattice with various values of d/X. 
The arrows indicate the cutoff magnetic fields associated with 
the vortex core. The respective cutoff values are 9.6, 7.9, and 
3.5 Bo. 
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FIG. 4: p(B_l/Bo) for triangular lattice for various values of 
d/X. The respective cutoff values are 9.2, 7.7, and 4.0 Bq. 



